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ABSTRACT 

Radiation pressure on dust grains may be an important physical mechanism driving galaxy-wide superwinds 
in rapidly star-forming galaxies. We calculate the combined dust and gas Eddington ratio (T) for the archetypal 
superwind of M82. By combining archival GALEX data, a standard dust model, Monte Carlo dust scattering 
calculations, and the Herschel map of the dust surface density distribution, the observed FUV/NUV surface 
brightness in the outflow constrains both the total UV luminosity escaping from the starburst along its minor 
axis (L*,uv) an d the flux-mean opacity, thus allowing a calculation of V. We find that L^ uv ~ 1 — 6 x 10 42 ergs 
s" 1 , ~ 2- 12 times greater than the UV luminosity observed from our line of sight. On a scale of 1-3 kpc above 
the plane of M82, we find that T ~ 0.01 -0.06. On smaller scales (~ 0.25-0.5 kpc), where the enclosed mass 
decreases, our calculation of L*,uv implies that T ~ 0.1 with factor of few uncertainties. Within the starburst 
itself, we estimate the single-scattering Eddington ratio to be of order unity. Thus, although radiation pressure 
is weak compared to gravity on kpc scales above the plane of M82, it may yet be important in launching the 
observed outflow. We discuss the primary uncertainties in our calculation, the sensitivity of T to the dust grain 
size distribution, and the time evolution of the wind following M82's recent starburst episodes. 
Subject headings: galaxies: formation, galaxies: general, galaxies: starburst, galaxies: individual: M82 



1. INTRODUCTION 

Galactic winds are crucial to the evolution of galaxies (see, 
e.g., Heckman et al. 1990). Several mechanisms have been 
proposed for driving them, including heating from supernovae 
(Chevalier & Clegg 1985) and stellar winds, heating from 
photoionization (Shapiro et al. 2004), and radiation pressure 
on dust grains (Murray et al. 2005, 2011, Zhang & Thomp- 
son 2012). Recent simulations indicate that radiation pressure 
may indeed be an important driver in galaxies with high star 
formation rates and gas densities (Hopkins et al. 2012). In 
the radiation pressure picture, radiation from hot stars drives 
dusty gas out of individual star-forming clouds, disrupting 
them, cutting off star formation, and then lofting this gas 
above the plane, where radiation from the entire galactic disk 
drives the gas fully out of the galaxy as a galactic superwind 
(Murray et al. 201 1, Hopkins et al. 2012). 

We examine the particular case of M82 (distance 3.63 Mpc, 
Gerke et al. 2011, Freedman et al. 1994), one of the nearest 
starburst galaxies to the Milky Way, in order to test whether 
radiation pressure is important in this case. There are exten- 
sive archival observations of M82 in many bands (see, e.g., 
McLeod et al. 1993, Westmoquette et al. 2007). However, 
the starburst core of M82, where most of the UV and opti- 
cal radiation that would be important for driving the wind is 
generated, is heavily obscured by dust along our line of sight 
because M82 is nearly edge-on (inclination: i = 80°, Mc Keith 
et al. 1993). The wind itself is quite bright in UV light on kpc 
scales above the plane (Figure [T}, and work by Hoopes et al. 
(2005) indicates that this emission is mostly due to dust in the 
wind scattering UV light from the central starburst into our 
line of sight, rather than other potential sources of in situ UV 
emission such as photoionization and shocked gas. 

The fundamental measure of whether radiation pressure is 
currently important in driving the wind is the Eddington ratio, 
T = Li/Ledd- To calculate V, we need to know the dynamical 
mass of M82 interior to our observation point a height z above 



the plane, M( < z), the flux-mean opacity of the dust/gas wind 
fluid k, and the true luminosity of the starburst that escapes to 
a height z from the disk, L*. 

Although the bolometric luminosity of M82 is fairly well 
known from FIR data (L IR = 5.9 x 1O 1O L , Sanders et al. 
2003), the shape of the SED in the UV above the plane, 
and specifically the UV luminosity L^.uv at a height z above 
the plane, is unknown without modeling. In principle, L*.uv 
could be as high as Lir. If so, assuming spherical symmetry, 

^Edd = 4irGcM * = 1.3 x 10" LoMio/ca 1 (1) 

where M w =M + /1O 1O M and k 3 = K/10 3 cm 2 g _1 of gas, 
normalized to a value appropriate for an SED dominated by 
UV light (Laor & Draine 1993). The Eddington ratio would 
then be r = L/L^dd ~ 0.5, suggesting that the radiation pres- 
sure force on the dusty gas could be comparable to gravity 
if the UV luminosity above the plane is of order Lir. For 
comparison, the total observed UV luminosity of M82 in the 
GALEX bands from Hoopes et al. (2005) is 1.2 x 10 8 L , or 
w 2 x 10 _3 Lir and « 9 x 10- 4 L Edd . 

By modeling the dust scattering, we work backwards using 
archival GALEX images of M82, presented by Hoopes et al. 
(2005) (see Figure [T}, to obtain the UV luminosity escaping 
perpendicular to the disk of M82. We then model the stellar 
population, which together with the scattered UV flux, allows 
us to obtain an SED of the radiation field at height z and to cal- 
culate the Eddington ratio with the chosen dust model. This 
calculation includes an MRN distribution (Mathis et al. 1977), 
STARBURST99 spectra (Leitherer et al. 1999), and dust grain 
properties calculated by Laor & Draine (1993). Changes in 
the grain size distribution are discussed in Section[3] 

In Section[2]we give an order of magnitude estimate of the 
UV luminosity that escapes perpendicular to the plane, as well 
as estimates of various parameters. In Section [3] we provide 
our results. In Section [4] we give a description of the uncer- 
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FIG. 1. — GALEX FUV (left) and NUV (right) images of M82. North is up and East is to the left. The red box on each image is our extraction aperture. 
Discussion of other extraction apertures and the implications for our results is provided in Section]?] M82 is at an inclination angle of 80° such that the southeast 
portion of the wind is tilted toward us. 



tainties in our work. A brief discussion and summary are pre- 
sented in Sections|5]and|6] 

While we were completing this work, Socrates & Sironi 
(2013) also presented a calculation of the Eddington ratio for 
M82 (top panel of their Fig. 1). However, their calculation 
makes use of the SED model of Silva et al. (1998), which 
makes no correction for the inclination dependence of the 
SED and the different shape of the radiation field emerging 
perpendicular to the disk. Our work is distinct in that we ac- 
count for the inclination dependence, estimate the SED shape 
from the scattered FUV and NUV light, and estimate k along 
the minor axis. In addition, by taking into account the ob- 
served rotation curve of the M82, we are able to estimate the 
Eddington ratio as a function of height along the minor axis. 
We are then able to make an estimate of the single-scattering 
Eddington limit within the starburst (Section[5]). 

2. ESTIMATES 
2.1. Order of Magnitude Calculation 

We expect that L+,uv < Lboi because at least a portion of the 
UV radiation in the starburst is absorbed within the system 
itself and never emerges to kpc scales above the plane. 

In order to estimate the UV flux which escapes the starburst 
region and scatters on grains a distance z above the plane, we 
make several simplifying assumptions: first, that the starburst 
disk is a point source; second, that the dust in the wind out- 
side of the starburst region is optically thin to the scattered UV 
(i.e., each photon scatters at most once); third, that the den- 
sity of the wind is constant along our line of sight at a given 
height above the plane of M82; and fourth, that the wind is 
illuminated as a right circular cone. Treating this system as 
steady-state, we can write down the radiative transfer equa- 
tion for this system: 



dl\ a\ L+\ 



ds 47T 47TT 2 



(2) 



where L+ \ is the specific starburst luminosity at wavelength 
A, r is the distance from the central source, p is the density 
of the dust, k s is its scattering opacity, and a is a wavelength- 
dependent factor that captures the anisotropy of dust scatter- 
ing in the UV. 



Because the total path length through the wind is compara- 
ble to the height above the disk, r changes significantly when 
integrating over a path through the wind that runs parallel to 
the disk. Along such a path (see Figure [2), 
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where s is the total thickness of the illuminated portion of the 
wind, s is the path length actually traversed by a given photon 
on this path, and z the height above the disk. This gives 
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The observed flux from a given aperture in the M82 wind is 
then 



I\ cos 6Wf2 = 2n I\ cos 9 sin 



(5) 



where the angle 9 C is the angular radius of the source aperture. 
For our work, we use an aperture with a radius of five GALEX 
pixels, or 7"5 (compare to GALEX's resolution of 4"2 and 
5 "3 in the FUV and NUV, respectively). Substituting, one 
finds an expression for the observed flux, which can then be 
inverted to give the luminosity of the starburst at wavelength 
A escaping the central starburst region 
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The GALEX images indicate that the dust is illuminated 
as a cone with a ~ 55° opening angle, so that s a ~ z. With 
reasonable values of the other parameters, the luminosity is 
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where Zk pc = z/kpc is the distance above the plane of 
M82, a 35 = a/0. 35 is a geometric factor that captures the 
physics of anisotropic UV scattering off dust grains, p_28 = 
p/10~ 28 gem" 3 is the mass density of the dust distribution, 
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FIG. 2. — Diagram of each of our model geometries for the dust distribu- 
tion. Blue shows the spherical model, and green the cylindrical model, r is 
the distance from the center of M82 to a given portion of the wind, s is the 
path length through the dust model to the observer from that portion, sq the 
total path length through the illuminated portion of the dust model, and z the 
distance from the midplane of M82, as defined in Section [2] The starburst 
core and extended optical disk of M82 (i = 80°) are also shown. Not to scale. 

^s.4.6 = ka/10 4,6 cm 2 g" 1 is the scattering opacity per gram 
of dust, Fobs,-n = -Fobs/10 -17 ergs s" 1 cm" 2 A is the observed 
flux, and A1516 = A/1516 A is the effective wavelength of the 
FUV band. 

The sum of the UV luminosity in the two GALEX bands 

is L^uv = AA fuv L*,fuv + AA NU vi*,Nuv ~ 2 x 10 42 ergss _1 , 
where AA is the effective bandwidth of the GALEX filters 
(268 A for FUV and 732 A for NUV). This is -5 times the 
total observed UV luminosity from Hoopes, et al. (2005) and 
about 10" 2 Lir. 

2.2. Parameter Estimates 

In order to estimate the density distribution of dust in the 
superwind, we use the dust column mass map from Figure 2 
of Roussel et al. (2010). This dust map is based on Herschel 
SPIRE 250 /im and 500 fjm data. Roussel et al. (2010) used 
these data to calculate a map of the dust mass surface density, 
E, with a resolution of 170 pc. They estimate that the total 
dust mass in the halo is ~ 1O 6 M and at least 65% of the dust 
is not associated with the superwind. From this map, we get a 
dust column density for each location in the wind. 

We assume two different model geometries: a spherical dis- 
tribution truncated at 5 kpc and a 5 kpc radius cylinder, both 
centered on M82 (Figures [2] and 13}. We then assume that the 
volume density of dust in each column is constant, so that 
p = £/s . In Figure[5]we show that both of these geometries 
produce dust density estimates within an order of magnitude 
of 10~ 28 gcm" 3 (Equation f?\. The overall shapes of the dif- 
ferent distributions are the result of the behavior of s a ; the 
spherical model produces higher densities further out because 
s„ rapidly shrinks at the extremities. It should be noted that the 
UV light from the starburst illuminates only a small portion of 
either of these models (FigurefTJ. This is consistent with the 
observed morphology of the UV and FIR maps, since the lat- 
ter shows an extended dust distribution, while the UV images 
show a smaller conical illuminated region. 

To calculate the dust opacity, we assume an MRN dust 




Distance from M82's Center (kpc) 

FIG. 3. — Dust density calculated from Roussel et al. (2010) for each of our 
model geometries along the extraction aperture shown in Figure[T]along the 
minor axis of M82 (positive and negative are north and south, respectively). 
The green curve is for the cylindrical model, and blue for the spherical (see 
Figure|2j. The central region is empty because the Roussel et al. (2010) map 
provides no data there (see their Figure 2). We applied a smoothing spline 
to the data for the figure. We also tested different extraction apertures, but 
found that using a different region of the wind had no significant effect on 
our results (Section]?}. 
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Distance from M82 nucleus (kpc) 

FIG. 4. — a as a function of z- The darker curves are for NUV, the brighter 
for FUV. The blue curves correspond to the spherical geometry, and the green 
the cylindrical. 

distribution (Mathis et al. 1977), with 50% each of graphite 
and silicate grains. We then use the absorption and scatter- 
ing cross sections from Laor & Draine (1993). We com- 
bine this with unreddened STARBURST99 models of the stel- 
lar spectra (Leitherer et al. 1999), and calculate a total flux- 
mean scattering opacity integrated over the entire SED of 
k.s.f ~4x 10 4 cm 2 g" 1 of dust (Equations |2j|7jl, assuming a bi- 
modal stellar population with components aged 4 and 9 Myr 
(Forster Schreiber et al. 2003). 

We use three dimensional Monte Carlo scattering simula- 
tions to estimate a. Each simulation run (one for each model 
geometry and bandpass) used 5 x 10 7 photons with a weight- 
ing scheme that guaranteed each photon scattered once (Wood 
et al. 2001). We assumed that the wind was optically thin, and 
that attenuation due to absorption was minimal and could be 
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FIG. 5. — Observed flux as a function of projected distance from M82 for the GALEX FUV (left) and NUV (right) data and our simulations. We used a circular 
aperture five GALEX pixels, or 7"5, in radius to calculate the average observed flux, F b s \ , at eachpoint. Negative values on the x-axis are for areas south of 
the plane of M82, while positive ones are to the north. The black curve is the GALEX data (FigurefTT. The dashed green and blue curves show results from our 
Monte Carlo simulations assuming the dust density distribution from Roussel et. al (2010), calculated assuming the cylindrical and spherical dust geometries, 
respectively, shown in Figure|2|and[3] Each of the simulation curves are normalized to the GALEX data at 2 kpc south of the plane. The southern portions of the 
curve are brighter due to the amsotropy in the scattering. The hole in the middle corresponds to the hole in the Roussel et al. (2010) dust density data; we simply 
set the density in that region to zero. The rapid falloff at the outer edges occurs for much the same reason; the simulation region abruptly cuts off at 4.8 kpc. 



ignored. To substantiate that assumption, we calculated the 
optical depth along each sight line. In the region for which we 
have dust density data, the total optical depth in either band 
for an average sight line is ~ 0.05 in the FUV, given the Rous- 
sel et al. (2010) dust density data and our model geometries. 
We ignored light from M82's outer disk, and approximated 
the starburst as a point-source. Modeling the starburst as a 
disk of radius 250 pc (Forster Schreiber et al. 2003) is a ~ 5% 
effect at 1 kpc from the plane, and becomes less important at 
larger distances. As shown in Figure[3] we have no dust den- 
sity data within ~ 1300 pc of the center of M82, because the 
Roussel et al. (2010) map provides no estimate of £ in that 
region due to difficulties in subtracting the disk PSF For the 
purposes of the simulation, we set the density to zero there, 
allowing photons to propagate freely outward until they en- 
counter the overlying wind material. This hole in the dust 
distribution results in zero scattering events into our line of 
sight. We ultimately removed the region from 1 .5 kpc south 
to 1.8 kpc north (see Figuresl5]and[8]). 

The geometric factor a is the ratio of the predicted flux with 
the Draine (2003a) phase function to the predicted flux for the 
case of isotropic scattering. Figure |4] plots a as a function of 
distance from the plane of M82. For the FUV band, a ~ 0.35 
and for the NUV band a ~ 0. 1 8 for both of our model geome- 
tries for the northern lobe of the wind. Due to the fact that the 
southern lobe is pointed towards us, a varies by a factor of ^5 
over that portion of the wind, from ^0.3 at 4.5 kpc to ~1.5 
at 1.5 kpc for the FUV band. This difference in the behavior 
of a between the northern and southern lobes arises from the 
shape of the Draine phase function. The Draine phase func- 
tion is relatively flat at the scattering angles covered by the 
northern lobe of the wind, but begins to rise more steeply as 
the dust moves into the forward-throwing regime. 

3. RESULTS 

3.1. The Intrinsic UV Luminosity 

We measured the observed flux from the wind along a lin- 
ear aperture perpendicular to the plane of M82 that extends 
5 kpc north and south of 9 h 55 m 52 s .7, 69°40'46" (Jackson et 



al. 2007). At each point along the line, we then measured the 
flux within a circular aperture five GALEX pixels in radius 
(GALEX's resolution is ~ 3 pixels). The average flux in each 
aperture and in each band was calculated and then corrected 
for Galactic extinction (E(B-V) = 0.159; Schlegel et al. 1998) 
in the GALEX bands (Wyder et al. 2007). The resulting flux 
in both bands is F a b s (see Equationj6l) and is shown as the solid 
lines in the two panels of Figure |5] 

The simulation results for the two adopted geometries are 
plotted against our extraction aperture through the GALEX 
data in Figure [5] Aside from an overall normalization applied 
such that all curves intersect at 2 kpc south, we have not tried 
to fit the GALEX observations. These curves use the two as- 
sumed wind geometries (Figure [3J together with the Roussel 
et al. (2010) dust data (Figure [3J. The overall difference in the 
flux between the northern and southern parts of the simulation 
results is due to the anisotropy in the scattering: the southern 
lobe of the wind is pointed towards us (i = 80°), and the dust 
is forward-throwing in the UV. 

Combining our mo del geometries, scattering calculations 
(for a(z); see Section |2~2] >, and the Roussel et al. (2010) dust 
data then gives L*,fuv and L*,nuv (Equation |6]), which we 
combine to give L+,uv, the escaping luminosity in the GALEX 
bands, at each z. We plot our inferred uv (z) for each model 
geometry in Figure [6] Given the assumptions of our model, 
£*,uv(z) should be constant as a function of height since only 
one total UV flux escapes the starburst region and there is 
minimal absorption perpendicular to the disk outside of the 
starburst region. The fact that L*,uv as calculated is not con- 
stant with z indicates that our model fails to capture the true 
geometry of the dust distribution. In principle, one could use 
the constraint that L*uvCz) = constant to determine (or con- 
strain) the true dust density distribution along the line of sight 
at every z. A future work could iterate between the scatter- 
ing calculations and the dust distribution using the joint con- 
straints of F b s , the Roussel et al. data, and uv(z) = constant 
to determine the true dust distribution. 

For our purposes, it is sufficient to take the range of L^ uv(z) 
from the spherical and cylindrical geometries as a measure of 
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FIG. 6. — The total UV luminosity in the GALEX bands escaping the star- 
burst region (L^ uv ) as determined from EquationPTlin both model geometries 
(same as Figure|5J, incorporating a(z.) from the Monte Carlo scattering cal- 
culations, p(z) from Figure]!] and the UV surface brightness from Figures^ 
and|5] If the geometry was correct and the dust density distribution perfectly 
determined, we would expect the inferred value of u V to be constant as a 
function of z- 

our systematic errors and we find that L+ t uy(z) ~ 1 — 6 x 10 42 
ergs s" 1 , a factor of ~ 2- 12 times larger than the total UV 
luminosity inferred by Hoopes et al. (2005) from our line of 
sight. 

3.2. The UV-NIR SED 

We next constrain the flux in the optical and NIR escap- 
ing perpendicular to the disk by comparing L*,uv with the 
true unreddened GALEX band luminosities expected from 
the work of Forster Schreiber et al. (2003), who modeled 
the stellar population in the central starburst region from NIR 
spectroscopy (Forster Schreiber et al. 2001). We produced 
STARBURST99 models for instantaneous bursts of star for- 
mation corresponding to a combination of a 4 and 9 Myr-old 
starbursts with parameters given by Forster Schreiber et al. 
(2003): L bol = 6.6 x 10 10 L Q and M+ = 6.1 x 10 8 M Q for the 
combined bursts at the current epoch. We weight the STAR- 
BURST99 spectrum by the GALEX response curves for each 
band and calculate the GALEX FUV and NUV band lumi- 
nosities of this fiducial, unextincted STARBURST99 model. 
As shown in Figure [JJ the unreddened UV luminosity is 
~ 2 x 10 44 ergs s" 1 . Our derived UV luminosities are much 
smaller (1 -6 x 10 42 ergs s" 1 ; Figure^, indicating that a sig- 
nificant fraction of the UV luminosity of the starburst is ab- 
sorbed on small scales. 

We model the absorption of the starburst UV flux using a 
simple foreground screen of dust. The natural logarithm of 
the ratio of our derived L*, FU v to the fiducial STARBURST99 
GALEX band luminosity then yields Tpuv- Since L*.fuv (or 
NU V, or combined UV) varies as a function of z (Fig. [7] Sec- 
tion 3.1 1, we derive a different optical depth at every z. We 
find that rpuv ranges from 1 .6 to 4, with the larger values cor- 
responding to the lower L*.uv in Figure [JJ and vice versa. 

To determine the optical and NIR luminosities emerging 
from the starburst region, we use the extinction law from 
Calzetti et al. (1994) and adjust it to provide our calculated 
r FUV at 1516 A, the effective wavelength of the FUV GALEX 
band. We then apply the Calzetti et al. law to our fiducial 
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FIG. 7. — The solid gray line shows the underlying assumed SED of the 
stellar population, as determined from the population synthesis modeling of 
Forster Schreiber et al. (2003). The black solid line shows an example of 
the extincted model used to calculate the Eddington ratio in Equation [8] and 
shown in Figurel8l This model gives L*,uv = 2 X 10 42 ergs s~', which corre- 
sponds to tfuv ^"j for a foreground screen. Given each of the derived values 
of L* uv(z) in Figure|6|we choose tfuv appropriately, and then use the SED 
model to calculate T(z). 

STARBURST99 SED, producing a reddened spectrum which 
we use for our Eddington ratio calculation. In Figure [JJ we 
show an example of how this reddened spectrum compares to 
the fiducial SED for one choice of rpuv- 

3.3. The Eddington Ratio 

The Eddington ratio for the combined dust and gas wind 
fluid, integrated over an MRN distribution, the redenned 
STARBURST99 spectral model, and using a thin, uniformly 
bright disk for the geometry of the starburst and a spherically 
symmetric mass distribution for the disk and galaxy, is 



Ar 2 



4cGM+m H (R 2 + r 2 ) 

I a^ 2 L^ x [Q, + Q s (l-g)]dad\ 



(8) 



where A is the normalization of the MRN distribution (cal- 
culated via the method in Laor & Draine 1993), is the 
dynamical mass of M82 interior to r, is the mass of a hy- 
drogen atom, r is the distance from the center of M82, R s is 
the radius of the starburst region, a is the radius of a dust grain, 
[3 is the slope of the grain size distribution, \ is the specific 
luminosity from Equation [6] Q d is the absorption efficiency, 
2s the scattering efficiency, and g = (cos 6) the average of the 
cosine of the scattering angle. It should again be noted that 
L*.a is calculated at every z above the plane in Figure ^(Fig- 
ure [Sjl, even though if the wind is optically thin, L+ x should 
be c onstant as a function of height above the plane (Section 
|3.1| ), We also calculate Tpuv for each z and use this in cor- 
recting for the extinction in M82's disk so that the SED varies 
as a function of distance from the disk (as illustrated in Fig- 
ure UJ. Although the variation in L^x as a function of height 
is an artifact of our assumed dust geometries, it allows us to 
understand the range of systematic error in our calculation of 

r(z). 

In order to determine M*(r), we use the work of Greco et 
al. (2012) (GMT12), who measured the rotation curve of M82 
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FIG. 8 . — Eddington ratio as a function of distance north of the plane of M82 as calculated from Equation|8] The green curve uses the cylindrical geometry, and 
the blue curve the spherical. All curves use the FUV data only; the curves generated from the NUV data are not significantly different. The gray curve shows the 
Eddington ratio calculated using the fiducial STARBURST99 SED with no extinction. The dashed black curve extends our calculation inward using a constant 
i*.uv ~ 2 X 10 42 ergs s for rpuv = 3. The black dashed and grey curves are relatively smooth because they keep L± \ constant with z. If we employ the full 
range of values implied for uv in Figure|6] the maximum deviation in the blue and green curves is a factor of six, from 0.01 - 0.06, and the maximum range in 
the dashed curve goes as high as ~ 0.2. The middle region is not shown for the colored curves due to the hole in the Roussel et al. (2010) dust map (see Figures[3] 
and[3), as explained in Section [2~2] See Section|5]for more discussion of this region. 

of dust grains, p is the dust mass density, and K m is per gram 
of dust. Since Q lv oc k, p , and L* > oc (a\K s ,\)~ l (see Equa- 
tion |6|), T\ oc k, p . \ Iol x n s . \. As the slope of the grain size 
distribution steepens, Ki-p,uv/Ks,uv grows; for an MRN distri- 
bution, /3 = -3.5 and K rp ,uv/K s .uv ~ 2.5, while for (3 = -4.5, 
K rp,uy/ K s,uv ~ 4, as shown in Figure [9] For shallower grain 
size distributions, the opacity ratio changes by relatively small 
factors until the size distribution gets close to flat (Figure |9j). 
Holding the dust-to-gas ratio fixed, changing the size distri- 
bution does not change the dust mass calculated by Roussel 
et al. (2010) by more than a few tens of percent. This is be- 
cause the FIR opacity is insensitive to the power law form of 
the grain size distribution at fixed dust-to-gas ratio, changing 
by only a few tens of percent from f3 = -5 to /? = -2. This 
means that over a reasonable range of slopes, the Eddington 
ratio does not change by more than a factor of three, too small 
to affect our overall result. 



on large scales. We assume circular orbits and a spherically 
symmetric mass distribution to determine MJf). These as- 
sumptions break down within ~ 0.5 kpc due to the influence 
of the bar, as detailed in GMT12 (see also Westmoquette et 
al. 2007, 2009, 2012). 

The solid colored lines in Figure [8] show our results for the 
Eddington ratio as a function of distance from the plane of 
M82 assuming a dust-to-gas ratio of f dg = 0.01. Factoring in 
our large uncertainties (Section [4j, V ~ 0.01-0.06 for our 
model geometries, indicating that radiation pressure does not 
presently drive the starburst wind on z >kpc scales. Note 
that if the dust and gas were not hydrodynamically coupled, 
the dust grains would be super-Eddington, with V — > r//<j g ~ 
1-3. We return to this issue in Section|5] 

The Eddington ratio T increases on smaller scales (z < 
1 kpc) because M + (r) decreases while L + . x likely stays 
roughly constant until z is comparable to the starburst's verti- 
cal scale height ~ 30 pc (WeiB et al. 2005). The dashed line 
in Figure [8] shows T(z) on small scales assuming the GMT12 
rotation curve and a single value for L*.uv = 2 x 10 42 ergs s -1 . 
Section[5]gives more discussion of T on small scales. 

As seen in Equation IS] the functional form of T depends 
heavily on the choice ofgrain size distribution and the optical 
properties of the grains. The radiation pressure cross section, 
Qrp = (2a + Gs(l -g), and the radiation pressure opacity, Krp, 
are related by 

n d Tra 2 Q lp = pKrp (9) 
for spherical grains of radius a, where «d is the number density 



4. ADDITIONAL UNCERTAINTIES 

There are two main sources of uncertainty in our estimate 
of the gravitational acceleration on the dust grains: the actual 
mass pulling on them, and the fact that the mass distribution 
is not perfectly spherically symmetric. The first is a relatively 
small uncertainty, as GMT 12 determined the total dynamical 
mass of M82 within 4 kpc from the center to within about 
10%. Treating M82 as a thin uniform disk with M = 10 10 M 
and radius 1 kpc instead of a sphere/point mass with the same 
mass reveals that the gravitational acceleration for the disk 
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Grain Size Distribution Slope (j3) 

FIG. 9. — Dependence of the Eddington ratio, T, on the slope of the grain 
size distribution, (3. The solid curve is for the FUV band, while the dashed 
curve is for the NUV band. For an MRN distribution as used in this paper, 
/3 = —3 . 5 . ft = is a completely flat grain size distribution. 

model is weaker by a factor of 1 / \[2 1 kpc above the plane 
from simple Newtonian gravity calculations. Thus while the 
fact that M82's potential well is not spherically symmetric can 
have a fairly large effect, and it implies a higher Eddington 
ratio, the effect is not large enough to affect our overall con- 
clusions, since our fiducial estimate of T is ~ 0.02. 

The radiation pressure force on the grains depends largely 
on two things: their scattering and absorption properties, and 
the flux incident on them. The biggest determinant of the first 
is our choice of grain size distribution. It is known that the 
MRN dust model used here does not reproduce the extinction 
curve towards, e.g., molecular clouds even in the Milky Way 
(see, for example, Weingartner & Draine 2001, Cardelli et al. 
1989), or the extinction present in starbursts (Calzetti et al. 
1994); M82 dust could have a different size distribution and 
composition as well. However, by the analysis in Section [3] 
this likely would not affect the answer by more than a factor 
of2-3. 

Our calculation of the flux incident on the grains is itself 
subject to several sources of uncertainty. As seen in Equa- 
tion^ it depends on the density of dust grains and their optical 
properties {L\ oc (a\pn s , a) -1 )- There are also three sources of 
reddening and extinction that go into calculating the flux: ma- 
terial in the disk of M82 and the starburst itself, extinction in 
the wind, and material in our own galaxy. We correct the last 
using the reddening corrections from Wyder et al. (2007), as 
described in Section [3] We assumed that the wind was opti- 
cally thin; we showedthat this was true for th e por tion of the 
wind we have dust density data for in Section [272] However, 
our results indicate that the starburst itself is optically thick 
in the UV (r FUV ~ 1.6-4). This large uncertainty in t fuv 
introduces a similarly large uncertainty in T. 

We correct for this reddening using the Calzetti et al. (1994) 
extinction law, as described in Section [3] This correction in- 
creases r by a factor of up to ~ 1.5. However, the Calzetti et 
al. extinction law assumes the dust is in a foreground screen, 
while the results of Forster Schreiber et al. (2001) indicate that 
modeling the dust and stars as intermixed is a better descrip- 
tion of the M82 starburst. In a series of comparisons between 
the shape of the emergent SED in the case of a foreground 
screen (e~ T ) and a mixed disribution of sources and absorbers 



([1 -e" r ]/r), we find that the Eddington ratio might change 
by a factor of ~ 1.5-2 at fixed L*uvfe). 

For the source geometry, we model the starburst as a point 
source in our order of magnitude calculations and our Monte 
Carlo simulations. Past ~1 kpc from the plane, this is a ^5% 
effect at most, so the uncertainty introduced in a\ is small 
compared to our other uncertainties. We also modeled the 
outer disk of M82 using STARBURST99 spectra as a 0.5 Gyr- 
old population of 4 x 1O 9 M (Mayya et al. 2006). Even as- 
suming no reddening or extinction between the starburst and 
the superwind, the outer disk in this model provides ~1% of 
the flux illuminating the wind, and therefore does not con- 
tribute significantly to the Eddington ratio. 

We also assume that the dust density is constant along the 
line of sight for a given z. While this cannot be strictly true, 
this assumption does not have a large effect on our calcula- 
tions. Although Equation [8|does contain p through L* \ and 
Equation[6j from EquationGlwe see that the observed flux of 
scattered light depends on the observed dust column, which is 
taken from Roussel et al. (2010). Thus, the main factor that 
will change given p(s) is a. We do not expect this change in 
a to heavily affect our results. 

In addition to the extraction aperture through the GALEX 
data we describe in this paper (Figure [TJ, we also examined 
two other extraction apertures through the wind, one centered 
a few tens of arcseconds northeast of our first one, and another 
centered a few tens of arcseconds southwest. We found that 
there was no difference in the overall conclusions; neither the 
dust densities nor the UV fluxes varied by more than a factor 
of two. Only the detailed shapes of the distributions changed. 

5. DISCUSSION AND CONCLUSIONS 

Our results indicate that even though the escaping UV lu- 
minosity is high compared to previous observations (Hoopes 
et al. 2005), the overall escape fraction of UV light is only 
~ 1 - 15%. Given that the dust opacity continues to rise to- 
wards the Lyman edge (Draine 2003b), it follows that there 
must be a low escape fraction of ionizing photons, likely less 
than ~ 10%. This is consistent with results obtained by sev- 
eral other groups (see, e.g., Heckman et al. 2001, Bergvall et 
al. 2006, Grimes et al. 2009, Heckman et al. 2011, Leitherer 
et al. 1995). Heckman et al. (2001) found that galactic winds 
do not necessarily clear paths through the ISM for ionizing 
radiation to escape. It does not appear that M82's superwind 
has cleared such a path. 

In this work, we generally treated the dust and gas as sepa- 
rate, except when we calculated the final Eddington ratio us- 
ing Equation [8] Ignoring magnetic fields and grain charging, 
which enhance dust-gas coupling, the scale over which grains 
of size a are hydrodynamically coupled to the gas is (Draine 
& Salpeter 1979, Murray et al. 2005) 

Am ^ lOao.i^nT'pc (10) 

where ao.i = a/(0.1 pm), n\ = n/cm" 3 is the gas number den- 
sity, and pj, = p/3g cm" 3 is the individual grain density. As 
shown by Figure [3] typical gas densities in the wind are 
0.1 - 0.01 cm" 3 , assuming a gas-to-dust ratio of 100, our 
model geometries for the dust density, and the Roussel et al. 
(2010) data. Thus, Am ranges from ~ 0. 1 — 1 kpc for 0.1 pm 
grains, and ~ 10-100pc for smaller grains. This means 
that the dust and gas are likely hydrodynamically coupled in 
the wind on large scales, and that the wind is sub-Eddington 
(T» 0.01-0.06 forz> 2 kpc). 
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Although our work indicates that the M82 superwind is 
sub-Eddington on large scales, there is still the question of 
whether it is super-Eddington on small scales (z ~ 0.25 - 
2kpc). If an MRN dust model and the rest of our assump- 
tions (see Section Hi are valid inside the region where we 
have dust density data from Roussel et al. (2010), then we 
can extend the Eddington ratio calculation inwards, given a 
determination of the enclosed dynamical mass. The rotation 
curve of GMT 12 implies M*(r < 250 pc) ~ 2-6 x 10 s M Q , 
but on these small scales their calculation is affected by the 
stellar bar in M82. Forster Schreiber et al. (2001) estimate 
M+{r < 250pc) ~ 8 ± 2 x 10 8 M from molecular and ionized 
gas tracers, similar, but somewhat higher than a naive appli- 
cation of GMT 12. Given these uncertainties, we estimate that 
for z ~ 250 pc 



r~o.i 



i*.UV 



8xl0 8 M Q \ 
2x 10 42 ergs s" 1 J VMJr < 250pc) ) 



(11) 



where there are factor of ^2-3 uncertainties in both L+ uv 
and M-t(r < 250pc). The dashed curve in Figure [8] shows an 
explicit calculation of Y on small scales assuming the GMT12 
mass estimate and L*,uv = 2 x 10 42 ergs s" 1 . Note that using 
the value of L*,uv — 6 x 10 42 ergs s" 1 (near the maximum in 
Figure [6]i increases the dashed curve in Figure [8] by a factor 
of 3 at equivalent tuv- Thus, this comparison indicates that 
the M82 superwind is currently sub-Eddington on z — 0.25- 
0.5 kpc scales. 

Another interesting possibility is that the superwind was 
super-Eddington in the past when the starburst was younger 
and brighter. The models of Forster Schreiber et al. (2003) in- 
dicate that the luminosity of the starburst was higher by about 
a factor of four roughly 6 Myr ago, for a bolometric luminos- 
ity of 2.6 x 10 11 L Q . Assuming a similar dust and gas dis- 
tribution to today's results in T ~ 0.05-0.3 on large scales 
and T ~ 0.2- 1 on smaller scales. However, if the gas surface 
density within the starburst was higher at that earlier time, this 
would lower the flux-mean opacity in the outflow, decreasing 
T throughout the wind. 

Finally, we consider the current Eddington limit for the 
dusty gas within the starburst itself. In this case, the appro- 
priate Eddington limit is the single-scattering limit discussed 
in Murray et al. (2005), Thompson et al. (2005), and Andrews 
& Thompson (201 1). If all of the radiation is absorbed or scat- 
tered once, then the single-scattering Eddington flux is ~ 
27rGS*SgC for a thin disk, where £* and Eg are the stellar 
and gas surface densities. Scaling to parameters appropriate 
to M82's nucleus, F Edd - 4 x 10 11 L Q kpc _2 (/ g /0.1)(M i /8 x 
10 8 M o ) 2 (250pc/fl) 4 [] where f g = £ f /E* is the gas frac- 
tion, whereas the observed FIR flux of M82 is Ffir — 3 x 
10 11 Lq kpc" 2 (L F iR/6 x 10 10 L©)(250pc/fl) 2 . Within the un- 
certainties, it remains plausible that M82's nuclear starbust re- 
gion currently exceeds the single-scattering Eddington limit, 
and may have significantly exceeded it 6 Myr ago, but clearly 
further work is required. Lastly, we note that many of the su- 
per star clusters in M82 exceed the single-scattering Edding- 
ton limit (Krumholz & Matzner 2009; Murray et al. 2010; 

1 Note that our use of f s = 0. 1 is prima facie inconsistent with tuv — 3 for 
the obscuring screen used to model the emergent SED in Section[3] since this 
UV optical depth would imply Eg ~ few X 1(T 3 cm 2 g -1 of gas. However, as 
discussed in Section [4] a mixed medium of sources and obscuration is more 
realistic, and implies tuv — 40 ± 30, which is consistent with f g = 0.1. 



Murray et al. 2011), and may thus be in part responsible for 
injecting gas into the hot outflow, or directly accelerating the 
gas to high velocities (Murray et al. 201 1; but, see KTumholz 
& Thompson 2013). 

6. SUMMARY 

We estimate the Eddington ratio for the dusty gas along the 
minor axis in the M82 superwind. We constrain the dust den- 
sity distribution using Herschel SPIRE data and the dust mass 
map provided by Roussel et al. (2010), and we use Monte 
Carlo simulations to verify the applicability of the MRN dust 
distribution and both of our model geometries in our cal- 
culations by comparing the simulation results with archival 
GALEX images of M82. We find that the combined UV lu- 
minosity in the GALEX bands that illuminates the superwind 
of M82 is 2- 12 times higher than that found by Hoopes et al. 
(2005) from our line of sight. 

We find that the dusty gas in the wind is highly sub- 
Eddington (T ~ 0.01-0.06) on vertical scales along the minor 
axis > 2 kpc (Figure 8J. On smaller scales, T is more un- 
certain, but increases because of the decrease in the enclosed 
mass. We estimate T ~ 0.05-0.3 at z ~ 250 kpc, with large 
uncertainties, depending on the dynamical mass of the central 
region, the distribution of absorbing material along the mi- 
nor axis, and the true UV luminosity that escapes the central 
starburst region (see Equation Hi. Note that if the dust were 



hydrodynamically uncoupled from the gas, which appears un- 
likely, the grains would be super-Eddington throughout the 
wind. In addition, the starburst was brighter by approximately 
a factor of four ~ 6 Myr ago (Forster Schreiber et al. 2003), 
which in principle could to lead to T as high as ~ 0.2 - 1 at 
z ~ 250 kpc. Finally, within the starburst itself, we find that 
M82 appears to be close to the single-scattering Eddington 
limit, implying that a significant fraction of the current ISM 
could be ejected by radiation pressure. 

Overall, our results provide quantitative constraints on the- 
oretical models of superwinds launched wholly or in part by 
radiation pressure on dust in dwarf starbursts like M82 (e.g., 
Murray et al. 201 1 ; Hopkins et al. 2012). 

There are a number of directions for future work. Perhaps 
most importantly, one could simultaneously model the scat- 
tering in a number of wavebands from the FUV through the 
NIR, using the full frequency-dependent dust scattering phase 
function. By iterating with the observed surface brightness 
profiles in each band, a self-consistent model might be de- 
veloped. Information could also be included on the MIR PAH 
surface brightness as a function of height. Such a study would 
constrain the dust geometry and grain size distribution more 
fully and perhaps allow for novel studies of the starburst pop- 
ulation in reflection off the wind. 
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